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The conventional use of two-body B decays to extract 7, although theo- 
retically clean, is currently statistics-limited. On the other hand, a bulk 
of data on three-body B decays is available from B factories. Applying 
the flavor-SU(3)-symmetric approach proposed in Ref. [2] to BaBar data, 
we find the highly promising result 7 = (8ll5(avg.) ±5(std. dev.))°. This 
establishes the use of three-body B decays as a viable alternative for the 
extraction of weak phases. In this preliminary analysis we have neglected 
several sources of uncertainties such as the effect of flavor-SU(3) breaking 
due to meson masses, and error correlations between input experimental 
parameters. A better understanding of these will improve the viability of 
this method. 

The Cabibbo-Kobayashi-Maskawa (CKM) phase 7 is conventionally extracted us- 
ing two-body B decays [Ij where the goal is to use the interference between the decay 
amplitudes for b — )• cus and h — )• ucs. This technique is theoretically clean because of 
the absence of penguins. However, current statistics on Cabibbo-suppressed modes 
of two-body B decays limits the precision of 7 obtained in two-body B decays. On 
the other hand, there is a fair amount of data available on B decays to charmless 
three-body final states from Belle and BaBar . However, extraction of weak phases 
from three-body decays is difficult because of two reasons. First, the charmless three- 
body final states do not involve distinct quark fiavors: one has to deal with penguin 
diagrams. Secondly, even fiavor-neutral three-body final states are not CP eigenstates 
: absence of indirect CP asymmetry makes it harder to extract the weak phase. It 
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was shown in Ref. [2] that in the hmit of flavor-SU(3) symmetry these difficulties 
can be overcome by combining Dahtz analyses of multiple three-body B decays and 
considering the part of the amplitudes that are fully symmetric under the exchange 
of any two final-state mesons. The purpose of this talk is to show an application of 
the method proposed in Ref. [2] using experimental data from BaBar . The theory 
details including the extraction technique are available in Ref. [2113] ■ 

We have used Dalitz analyses for the following processes available from BaBar : 
B° ^ KsTT+TT- P], KsK+K- PS, K+ir-TT^ [6], and 3Ks [7]. In addition it is neces- 
sary to include data from 2?+ — )• K+tt+tt" to estimate SU(3) breaking as described 
in Ref. [2]. However, our preliminary goal is to establish the viability of the method 
and we have ignored effects of SU(3) breaking. In order to relate the amplitudes of 
the different three-body decay modes they are written in terms of flavor-SU(3) dia- 
grammatics. In three-body B decays, neglecting power-suppressed contributions from 
weak-annihilation topologies, there are ten distinct flavor-SU(3) topologies, "Color- 
favored Tree" (Ti,T2), "Color-suppressed Tree" (Ci,C2), "Color-favored electroweak 
Penguin" {Pewi: Pewt), "Color-suppressed electroweak Penguin" {Pewi^ Pewi) 
"Gluonic Penguins" {Ptc,Puc)- 

The amplitudes for the relevant B decay processes, under the assumption of flavor- 
SU(3) symmetry may be expressed in terms of these twelve flavor-topology parame- 
ters. In order to construct a CP eigenstate from the three-body flnal state, Ref. [2] 
considers the part of the amplitude that is totally symmetric under the interchange 
of any two flnal state mesons. An added advantage of using the totally-symmetric 
flnal-state amplitude is that the "electroweak Penguin" topologies can now be related 
to the "Tree" topologies [3]. 

Although it is possible to measure both the direct and indirect CP asymmetries in 
B^ — 7- 3Ks, these observables are currently limited by the available statistics This 
can be implemented by setting P^c to zero. The three-body amplitudes may now be 
written in terms of "effective diagrams" (linear combinations of the flavor-topology 
amplitudes) as follows: 

A{B^ K^K°Tf),y^ = a , 
V2A{B° K+K^K-),y^ = -ce'^ - a + nb , 

v^A(5° ^ K\+7v-)sy„, = -de''' -a + Kd , (1) 
where in terms of the flavor-topology amplitudes a, b, c and d are: 

a = -Pte + «: (^^Ti + + iCa) , 
b = T, + C2 , c = T2 + Ci , d = T^ + Ci , . (2) 

The fully-symmetric three-body amplitudes are dependent on the kinematics of 
the final-state particles. This dependence can be conveniently represented on a Dalitz 
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plot. Dalitz analyses for the relevant three-body decays have been performed by 
BaBar . A convenient model used to represent the dynamics of the three-body final 
state is the Isobar model where the final state amplitude is written as follows: 

^DP = ATdp J2 Cj-e''^ (si2, S13) , (3) 
j 

where the index j represents an intermediate state, A/dp serves the purpose of nor- 
malizing the amplitudes to the measured rates for the relevant processes, Sik is the 
invariant mass of two final state mesons labeled i and k respectively, and Cj and 6j 
represent the amplitude and phase of the jth isobar coefficient. Although there are 
three choices for Sik, the sum of all three is a constant and hence only two of them 
are independent. 

The values of cj and 9j were taken from the respective BaBar publications. Al- 
though, the error bars on these quantities are correlated, such correlations are often 
not available in the literature. In our preliminary analysis we have chosen to ignore 
correlations, noting that a more complete work would necessarily include them. We 
are now able to compute the totally symmetric parts of the amplitudes as follows: 

Aym = -^iA{Si2, S13) + ^(Si3, S12) + A{Si2, S23) 

+A{S23, S12) + A{S23, S13) + A{Si3, S23)) . (4) 

A similar exercise can be done to obtain the symmetric part of the amplitude for 
CP-conjugate decay process. These may then be used to construct the observable 
quantities for every point on a given Dalitz plot as follows: 
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A simple count of parameters and observables tells us that one can extract 7 in 
this scenario. For every point on the Dalitz plot we have eight independent unknown 
parameters : the magnitudes and relative phases of a, b, c and d (one overall phase is 
arbitrary) and 7. However, we have nine observables : X^p from all four processes, 
Ydp from all but B — )■ 3Ks (due to low statistics CP-asymmetry measurements aren't 
available for this decay), Zdp from — )■ Ks{'n'~^'n'~ , K^K~). Note that the K^-k'ti^ 
final state is not fiavor neutral, and doesn't give us a Zpip. The advantage of this 
method is that we are able to determine 7 independently for every point on the Dalitz 
plot. Therefore an average over Dalitz plot points is expected to significantly reduce 
the error in 7. 

In Fig.[T]we show the kinematic boundaries of the Dalitz plots for the four relevant 
decay processes. In addition we show the relevant axes of symmetry, which divide 
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each Dalitz plot into six zones, ^sym is totally symmetric under the exchange of any 
two Sjfc's. Therefore to avoid overcounting we are restricted to one of the six zones. 
Within the chosen sixth of the Dalitz plots we pick fourteen points to use for the 
determination of 7, which have been displayed in Fig. [TJ We have chosen the points 
to lie on a grid with the squared momenta separated by 2GeV^ on either directions. 
This choice is completely arbitrary. For each point we perform a fit. 
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Figure 1: Dalitz boundaries and symmetry axes for the four relevant B decays. The 
blue dots represent the points that were used for extraction of 7, while the red dots 
represent the points close to the boundary of the Dalitz plots that were discarded. 

In Fig. [21 we plot the minimum as a function of 7 between and 180°, for two 
different points on the Dalitz plots. The first plot shows one solution with 7 ~ 90° 
which can be clearly distinguished from the other solutions at the Ax^ = 1 level, and 
several other indistinguishable solutions closer to 7 ~ 30°. The second plot corre- 
sponds to a point close to the Dalitz boundaries and shows multiple indistinguishable 

minima. We find similar results for other points close to the Dalitz boundaries. 
These solutions require a better understanding of the underlying dynamics and have 
currently been left out while computing the final result for 7. We find that a clearly 
distinguishable solution exists for nine points that are further away from the Dalitz 
boundaries. 

In Tabled] we present our results for 7 obtained from each Dalitz plot point where 
we found a distinguishable solution. A simple mean over the values of 7 quoted in 
Table [U gives us the following result: 

7 = (81+^(avg.) ± 5(std. dev.))° , (6) 
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Figure 2: minimum as a function of 7 for two different Dalitz plot points. The 
plot in the left panel shows one clearly distinguishable solution for 7 ~ 80°. The plot 
on right shows multiple indistinguishable solutions. 



where the first error was obtained by square averaging over the individual errors and 
the second error denotes the standard deviation from the mean of the central values 
of 7 for the nine Dalitz plot points used. 

The error bar on 7 obtained above is significantly smaller than that currently 
obtainable in two-body decays. However, we have neglected a few sources of errors. 
For example we have ignored correlations between isobar coefficients. In practice such 
correlations may introduce additional error. Furthermore we have assumed flavor- 
SU(3) symmetry. However, the masses of the daughter pions and kaons as well as 
the intermediate resonances break fiavor-SU(3) symmetry. A direct consequence of 
fiavor-SU(3) breaking is that the Dalitz boundaries for a final state K-rni is different 
from that for KKK. As a consequence points very close to the Dalitz boundaries 
don't yield distinguishable solutions for 7. One way of alleviating this problem is by 
including many more points on the Dalitz plot for evaluating 7. Note that a more 
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(4, 4) 
(6, 2) 
(6, 4) 
(6, 6) 
(8, 2) 
(8, 4) 
(8, 6) 


70+12 
80li° 

80^}^ 
781JI 


(8, 8) 
(10, 2) 
(10, 4) 
(10, 6) 
(10, 8) 
(12, 2) 
(12, 4) 


80±1^ 

861}! 
791}^ 
8lli° 

8811^ 



Table 1: 7 extracted from several points within the Dalitz plots. 
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complete analysis of the combined data set for the four Dalitz plots can avoid errors 
due to the neglect of correlations between isobar coefficients. In order to get an idea 
of how large the effect of SU(3) breaking is, it would be interesting to include effects 
of SU(3) breaking due to one extra parameter. 

The work presented in this talk was aimed at motivating interest toward a new 
method for extracting 7 from three-body B decays. In this preliminary analysis we 
have neglected several sources of errors, however, the value of 7 obtained is quite 
promising. The hope is to find an alternative viable method for extraction of 7 using 
the already available statistics on three-body B decays from B-factories. 

This work was financially supported by the NSERC of Canada (BB, DL) and 
FQRNT of Quebec (Ml). BB would like to thank the conference organizers for a 
wonderful stay in Cincinnati. 
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